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SCALING  THE  ENERGY  SPECTRA  OF  UNDERWATER  EXPLOSION  SHOCK  WAVES 

By 

Ermine  A.  Christian 


ABSTRACT:  Scaling  factors  with  which  the  energy  spectrum  for 
an  underwater  shock  wave  can  be  scaled  from  one  set  of  conditions 
(i.e,,  explosive  composition,  charge  weight,  and  range)  to  a 
second  set  of  conditions  are  derived.  Analogous  scaling  factors 
for  an  ideal  acoustic  wave  are  also  shown.  The  spectrum  measured 
at  1  yard  from  a  1-lb  charge  is  scaled  to  two  other  conditions 
(viz.,  100  yards  from  a  1-lb  charge,  and  280  yards  from  a  1.8- 
lb  charge)  for  which  measured  values  are  also  available.  In 
both  cases  the  explosion  scaling  is  in  agreement  with  the  data 
whereas  the  ideal  acoustic  scaling  gives  spectrum  levels  that 
are  too  high  at  the  higher  frequencies. 

The  allowable  ranges  of  variables  for  the  spectrum  scaling 
functions  shown  here  are  estimated  to  be: 

0.5  <  frequency  (kc)  <  20 

4  x  lo”1*  s  W^/^/R  s  3,  where  W  is  charge  weight  in 
pounds  and  R  is  range  in  yards. 

Details  of  the  spectral  analysis  of  measured  shock  waves 
and  a  method  of  estimating  the  effect  of  explosive  composition 
on  spectrum  level  are  given  in  appendices. 
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shock  wave  energy  spectra  for  different  sets  of  experimental 
variables  or  to  scale  spectra  from  one  set  of  conditions  to 
another.  In  this  report  generalized  scaling  factors  for 
exponential  explosion  pulses  are  shown,  end  the  differences 
between  the  scaling  laws  for  explosion  shock  wave  and  ideal 
acoustic  pressure  waves  are  discussed. 
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SCALING  THE  ENERGY  SPECTRA  OP  UNDERWATER  EXPLOSION  SHOCK  WAVES* 


1 .  INTRODUCTION 

1.1  At  very  short  ranges  from  an  underwater  explosion  the 
relationships  between  charge  weight,  range,  and  explosive  compo¬ 
sition  have  been  well  established  from  numerous  measurements  of 
the  damage  potential  of  a  charge  fired  underwater  (ref.  (1))**. 

A  number  of  measurements  have  also  been  made  at  great  distances 
from  small  charges  used  as  acoustic  sources  (e.g.,  refs.  (2),  (3), 
and  (9)),  and  it  has  been  found  that  at  ranges  of  many  miles  8- 
coustic  wave  propagation  laws  adequately  describe  the  shock  wave 
propagation.  At  Intermediate  ranges,  however,  within  a  few  miles 
of  a  charge  weighing  several  pounds,  the  relationships  between 
explosion  and  acoustic  waves  have  not  been  well  defined.  Theory 
predicts  that  the  shock  wave  propagation  laws  will  approach  asymp¬ 
totically  the  ideal  acoustic  laws,  but  Arons  (ref.  (K\)  found 
that  the  near-field  similitude  laws  described  the  explosion  wave 
form  out  to  greater  ranges  than  expected.  It  is  this  intermediate 
region,  which  is  "short-range"  from  the  point  of  view  of  most 
previous  explosion  acoustics  studies  and  is  "long-range"  from 

the  point  of  view  of  explosion  damage  studies,  that  is  of  primary 
interest  for  explosives  echo  ranging  application. 

1.2  Within  thlB  intermediate  region  acoustic  approximations 
are  suitable  for  some  of  the  shock  wave  characteristics  but  are 
significantly  in  error  for  others.  For  example,  the  propagation 
velocity  of  the  wave  front  decreases  from  its  initial  high  value 
to  essentially  acoustic  velocity  within  a  distance  of  perhaps  20 
charge  radii  away  from  the  detonation.  The  shock  wave  energy 
spectrum,  on  the  other  hand,  reflects  the  non-linearities  of  the 
initial  high-amplitude  wave  out  to  ranges  far  beyond  the  point 
at  which  the  wave  velocity  is  acoustic. 

1.3  The  spectrum  of  the  shock  wave  is,  of  course,  a  function 
of  the  charge  weight.  It  is  also  a  function  of  range,  and  for 
discrete  frequency  components  of  the  broad  band  pulse  the  range 
dependency  is  different  from  the  range  dependency  of  a  sinuoidal 
pulse  of  the  same  frequency.  Because  the  peak  pressure  falls  off 
with  range  at  a  greater  rate  than  1/R  (the  ideal  scoustic  wave 
decay)  and  the  profile  of  the  wave  continues  to  broaden,  the 
spectral  distribution  of  the  shock  wave  changes  as  the  pulse 
propagates  outward.  These  effects,  which  are  not  present  In  Ideal 
acoustic  waves,  continue  out  to  surprisingly  great  ranges. 

resume'  of  this  paper  was  presented  at  the  Sixty-Fourth  Meet¬ 
ing  of  the  Acoustical  Society  of  America,  November  1962,  Seattle, 
Washington, 

♦♦References  are  listed  at  the  end  of  the  report. 
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1.4  Since  detailed  data  on  shock  waves  several  miles  from 
small  charges  are  sparse,  it  ia  frequently  necessary  to  scale  an 
available  spectrum  to  other  conditions  of  interest.  It  is  clear 
from  the  above  that  if  the  spectrum  is  scaled  as  though  the 
shock  wave  were  propagating  as  an  ideal  acoustic  wave,  errors  in 
scaled  spectrum  levels  can  result.  The  magnitudes  of  such  errors 
depend  upon  the  charge  weights,  the  distances,  and  the  frequencies 
Involved.  The  relationships  between  explosion  wave  and  ideal 
aooustio  wave  treatments  of  the  shock  wave  energy  spectrum  are 
the  subject  of  this  paper. 

1.5  In  the  so-called  "explosion  wave"  treatment  here  the 
shock  wave  energy  spectrum  level  as  a  function  of  charge  weight 
and  range  ia  determined  in  the  manner  of  Weston  (ref.  (5)).  The 
shock  wave  is  assumed  to  be  a  steep-fronted  exponential  wave, 
the  peak  pressure  and  decay  constant  of  which  are  given  by  ex¬ 
plosion  similitude  laws  that  describe  the  shook  wave  near  the 
charge.  The  "ideal  acoustic  wave"  treatment  assumes  that  the 
pressure  falls  off  with  distance  as  1/R  and  that  there  is  no 
change  in  the  pulse  shape  as  the  wave  propagates.  For  both 
treatments  shock  wave  energy  spectrum  functions  are  presented 

in  the  form  of  scaling  relationships  with  which  the  spectrum 
for  one  set  of  experimental  conditions  (i.e.,  explosive  compo¬ 
sition,  charge  weight,  range)  can  be  scaled  to  a  second  set  of 
experimental  conditions. 

1.6  Although  only  the  acoustic  wave  treatment  is  denoted 

aa  "ideal",  the  shook  wave  treatment  does,  in  fact,  imply  several 
idealizing  assumptions.  The  range  out  to  whioh  similitude 
relationships  obtain  la  open  to  question,  since  measurements  of 
very  low  amplitude  waves  are  highly  subject  to  instrumental 
errors.  In  addition,  the  use  of  an  exponential  wave  form  implies 
the  assumptions  that: 

(1)  The  rounding  of  the  front  due  to  dissipative  effects 
(ref.  (6))  does  not  decrease  the  spectrum  level, 

(2)  The  contribution  of  the  slowly-decaying  pressure 
in  the  tall  of  the  shook  wave  ia  negligible,  and 

(3)  The  wave  is  propagating  in  an  infinite,  homo¬ 
geneous  medium  ao  that  the  pulse  shape,  and  hence  the  spectrum 
level,  ia  not  affected  by  either  reflections  (ref.  (7))  or 
refraction  (ref,  (8)), 

A  comment  concerning  the  range  of  variables  over  which  these 
assumptions  are  reasonable  la  included. 
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2.  SHOCK  WAVE  ENERGY  SPECTRUM  SCALING 

2.1  For  an  exponentially  decaying  pulse  the  Fourier  Integral 
for  the  energy  spectrum  level  Is: 

E(f )  «  ~  [  ]  (1) 

»° L  (i«V) J  w 

where  p  is  the  Initial  pressure  of  the  wave,  e  Is  the  decay 
constant,  pc  is  the  aooustlc  Impedance,  and  u  is  the  angular 
frequency,  cu  «=  2 ir  times  the  frequency  In  cps.  For  an  explosion 
shock  wave,  the  pressure  and  time  eonstfcnt  In  Equation  (l)  can 
be  determined  from  similitude  equations  of  the  form: 


p  -  k  (wV3A)a 

(2) 

6  -  m  (W^/R^.W1/3 

(3) 

Here  W  is  charge  weight,  R  Is  distance  from  the  detonation,  and 
the  coefficients  and  exponents  are  functions  of  the  charge  compo¬ 
sition.  Combination  of  Equations  (1),  (2),  and  (3)  leads  to 
the  following  general  scaling  factors  that  can  be  used  to  scale 
the  shock  wave  energy  frequency  spectrum  E^ff^)  for  one  Bet  of 

conditions  to  the  speotrum  Eg(fg)  for  a  second  set  of  conditions: 
Eg  -  K£E1  when  fg  -  Kff1,  where 


2,2  For  a  given  explosive  material  the  above  spectrum 
scaling  factors  reduce  to  simple  forms  for  two  cases  that  are 
frequently  of  interest: 
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(a)  Constant  range,  varying  charge  weight 


2<#  +  8  +  ^  Kf(w)  .  $  ]‘(s  +  1)/3 


(6) 


(b)  Constant  charge  weights  varying  range 


-[^]2(a+ B)i  -  ty’ 


When  both  charge  weight  and  range  are  varied,  the  products  of 
these  factors  are  used, 

2,3  Since  it  has  been  found  that  the  values  of  a  0  1.13  and 
0  «  -0.22  determined  experimentally  for  TNT  (ref.  (4;)  are  typioal 
of  the  exponents  for  most  high  explosive  materials*,  these  values 
are  used  here  to  compare  explosion  wave  and  Ideal  acoustic  wave 
scaling  factors.  For  an  ideal  acoustic  wave  there  is  no  spread¬ 
ing  of  the  profile  (l.e.,  0  -  0)  and  the  pressure  decays  as  1/R 
(l.e.,  o  ■  1.0) .  In  Table  X  the  spectrum  soallng  factors  of 
Equations  (6)  and  (7)  for  TNT  and  for  an  ideal  acoustic  wave  are 
summarized. 


TABLE  1,  SPECTRUM  SCALING  FACTORS 


Constant  Range,  Constant  Weight, 

Varying  Weight _ Varying  Range 


TNT 

W9  1,27 

W5  -0.26 

R,  1.82  R.  0.22 

(a  -  1.13,$  -  -0.22) 

Kg  -  (y-) 

iKf 

“e  iKf 

Ideal  Acoustio 

1.33 

W2-0.33 

Rn  2.0 

(«  -  1.0,  0  -  0) 

;Kf  -  (y-) 

%  *  <HJ>  -  1.0 

The  ideal  acoustic  factors  of  Table  1  correspond  to  the  spectrum 
scaling  constants  given  by  Weston  in  reference  5.  In  figures 
1  and  2  the  functions  of  Table  1  are  shown  graphically.  Dashed 
lines  in  these  figures  represent  the  TNT  scaling  factors  and 
solid  lines  the  ideal  acoustic  wave  factors.  For  convenience, 
the  values  of  Kg  are  shown  in  decibels  (10  log  Kg)  so  that  they 

**The  effect  of  charge  composition  on  spectrum  level  is  discussed 
In  Appendix  2. 
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can  be  simply  added  to  or  subtracted  from  a  known  spectrum  level 
that  Is  in  dB. 

2.4.  In  figure  1  the  scaling  factors  of  Equation  (6),  for  the 
case  of  a  oonstant  range  and  varying  charge  weight,  are  shown. 
Here  the  curves  are  normalized  to  a  charge  weight  of  one  pound, 
so  that  If  the  spectrum  for  a  one-pound  charge  Is  known  at  some 
particular  range,  the  spectrum  for  other  charge  weights  at  that 
same  range  can  be  determined  from  soale  factors  read  directly 
from  figure  1.  For  example,  using  explosion  factors  the  spectrum 
(E2(f2))  of  a  10-pound  charge  (at  the  same  range)  Is  obtained 
from  the  spectrum  of  a  one-pound  charge  (E1(f1))  by  adding  12.7 

dB  to  the  energy  levels  and  multiplying  the  frequencies  by  a 
factor  of  0.55*  In  other  words,  if,  for  a  one-pound  charge, 

E1  »  y  (dB)  at  f1  »  x  (cps),  at  that  same  range,  for  a  10-pound 

charge 


Eg  «  y  (dB)  +  12.7  at  f2  **  0.55  x  (cps). 

Similarly,  if  the  one-pound  spectrum  1b  being  scaled  to  that  of 
a  charge  weighing  only  0.01  pounds,  the  explosion  soaling  factors 
would  be: 


Eg  *»  y  (dB)  -  25.5  dB  at  f2  -  3.3  x  (cps) 

Thus,  at  any  given  range.  Increasing  the  charge  weight  Increases 
the  spectrum  level  and  shifts  the  spectrum  towards  lower  frequen¬ 
cies,  while  decreasing  the  oharge  weight  lowers  the  level  and 
shifts  the  speotrum  to  higher  frequencies.  As  can  be  seen  from 
the  two  sets  of  curves  In  figure  1,  the  differences  between  ex¬ 
plosion  and  ideal  acoustic  scaling  of  the  spectrum  arise  from  the 
frequency  soaling  factor,  rather  than  the  level  scaling  fa'ctor*. 
It  Is  also  apparent  from  figure  1  that  the  two  sets  of  scaling 
faotors  give  essentially  the  same  results  for  charge  weights 
between  about  0.1  and  100  pounds.  For  oharge  weights  less  than 
0.1  pounds,  however,  significant  errors  In  the  spectrum  could 
be  introduced  by  use  of  the  ideal  acoustic  factors;  the  magnitude 
of  the  error  would  depend  upon  the  slope  of  the  level  at  the 
particular  frequencies  of  Interest. 

2.5  In  figure  2  the  scaling  factors  of  Equation  (7)  for  the 
case  of  constant  charge  weight  and  varying  range  are  shown.  Here 
the  values  are  normalized  to  a  100-yard  range;  scaling  factors 
are  found  as  in  the  preceding  example.  Figure  2  shows  that  for 
a  given  charge  weight  the  spectrum  level  Is  decreased  and  shifted 

•The  greater  significance  of  the  frequency  factor  arises,  of 
course,  from  the  use  of  decibels,  rather  than  ratios,  to  express 
spectrum  levels.  Numerically  the  two  factors  are  about  equally 
significant. 
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towards  lower  frequencies  as  range  is  increased.  Thus,  although 
the  level  Is  lowered  by  either  a  decrease  In  weight  or  an  In¬ 
crease  in  range,  the  frequency  shifts  are  In  the  opposite 
directions  for  the  two  cases.  As  was  the  case  for  constant  range, 
the  difference  In  explosion  and  Ideal  acouBtic  scaling  for 
spectra  from  different  charge  weights  is  due  to  the  frequency 
shift  rather  than  the  level  change.  From  figure  2  it  is  also 
apparent  that  there  are  only  relatively  minor  differences  in  the 
two  sets  of  scaling  factors  for  ranges  between  about  10  and  1000 
yards.  However,  In  scaling  the  spectrum  at  100  yards  back  to 
that  at  1  yard,  the  reference  range  customarily  used  for  sonar, 
differences  between  the  two  methods  of  scaling  may  be  significant*. 
As  noted  previously,  the  magnitude  of  such  scaling  errors  varies 
with  frequency. 

2.6  A  final  example  of  the  use  of  figures  1  and  2  Is  given 
for  the  case  of  changing  charge  weight  and  range  simultaneously. 
The  spectrum  for  a  one-pound  charge  at  100  yards  can  be  scaled 
to  that  of  a  0.1-pound  charge  at  10  yards  with  the  following 
factors: 


10  log  Kg  due  to  charge  weight  decrease  is  -12.7  dB(Fig.l) 

10  log  due  to  range  decrease  is  +18.2  dB(Fig,2) 

so  that  the  level  is  increased  by  +'  £.£>  clB 

Kf  due  to  charge  weight  decrease  is  1.82  (Fig.  1) 

Kf  due  to  range  decrease  is  1.65  (Fig.  2) 

so  that  frequencies  are  multiplied  by  (1.82  x  I.65  »  3.0) 

2.7  The  curves  of  figures  1  and  2  are  directly  applicable  to 
weight  and  range  ratios,  as  well  as  to  the  values  of  weights 

and  ranges  shown.  Inother  words,  the  spectrum  of  a  5.76-pound 

charge  can  be  scaled  to  the  spectrum  of  a  0.0576-pound  charge 
with  factors  read  at  10“2  on  the  abscissa  of  figure  1.  Simi¬ 
larly,  factors  for  scaling  a  spectrum  from  50  yards  to  5000  yards 
can  be  read  at  101*  (since  the  curves  shown  are  normalized  to  10*) 
on  the  abscissa  of  figure  2.  It  should  also  be  noted  that 
the  curves  of  figures  1  and  2  can  be  used  to  scale  spectra 
between  any  two  values  of  weight  and  range,  even  though  the 

•Weston  thei\(5))  considers  100  yards  to  be  a  suitable  reference 
range  for  determining  the  acoustic  source  level  of  an  explosive 
charge,  since  finite  amplitude  effects  are  comparatively  small 
beyond  about  100  yards j  this  distance  Is  also  typical  for  the 
range  from  oharge  to  first  boundary  reflection.  There  is,  how¬ 
ever,  a  tendency  to  scale  Weston’s  100-yard  data  back  to  a  l-y8rd 
range,  e.g.,  reference  9. 
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appropriate  factors  cannot  be  read  directly  from  the  ordinate 
scales  shown.  For  such  use  the  scaling  values  that  are  needed 
are  those  which  would  shift  the  correction  factor  curves  shown 
here  so  that  10  log  Kg  -  0  and  Kf  -  1,0  for  the  weight  and  range 

of  the  known  spectrum.  For  example,  the  factors  for  scaling 
the  spectrum  at  one  yard  from  a  given  charge  weight  to  the 
spectrum  for  that  same  charge  at  1000  yards  can  be  obtained  from 
the  values  shown  for  these  two  ranges  in  figure  2.  In  this  case, 
using  explosion  factors,  spectrum  levels  would  be  decreased  by 
55  dB  (the  difference  between  the  +37  dB  value  at  1  yard  and  the 
-l8  dB  value  at  1000  yards)  and  the  corresponding  frequencies 
would  be  multiplied  by(0, 60/2,75  =>  0.22),  the  ratio  of  the 
factors  at  the  two  ranges. 

3.  COMPARISON  WITH  EXPERIMENTAL  RESULTS 

3.1  In  figure  3  the  energy  flux  spectrum  levels  at  one  yard 
and  at  100  yards  from  a  one-pound  TNT  charge  are  shown.  The 
uppermost,  solid  curve  was  obtained  by  digital  analysis  of  a 
shock  wave  recorded  2.5  ft  from  the  charge  and  scaled  to  a  one- 
yard  range  using  Equation  (7).  Experimental  conditions  and 
details  of  the  spectral  analysis  for  this  curve  and  the  measured 
spectrum  of  figure  4  are  given  in  Appendix  1.  The  100-yard  data 
(open  circles)  of  figure  3  are  the  octave  band  measurements  given 
by  Weston  in  Table  2  of  reference  5*.  The  two  dashed  curves  for 
the  spectrum  at  100  yards  were  scaled  from  the  1-yard  spectrum 
with  the  factors  of  figure  2;  the  long  dash  represents  the  ex¬ 
plosion  faotors  and  the  short  dash  represents  ideal  acoustic 
factors. 


3.2  The  measured  and  scaled  spectra  at  a  range  of  280  yards 
from  a  1.8-pound  charge  are  shown  in  figure  4,  The  solid  curve 
is  the  spectrum  computed  from  a  measured  shock  wave  and  the 
dashed  curves  were  scaled  from  the  1-pound,  1-yard  curve  of 
figure  3,  Again,  the  long  dash  Indicates  scaling  with  explosion 
factors  and  the  short  dash  Indicates  scaling  with  ideal  acoustic 
factors.  Since  both  charge  weight  and  range  were  changed  here, 
both  figures  1  and  2  were  used  for  the  scaled  curves  of  figure  4, 

3.3  At  both  100  yards  and  280  yards  the  explosion  wave 
scaling  factors  give  excellent  agreement  between  measured  and 
scaled  shock  wave  spectra;  except  at  the  highest  frequency  for 
the  100-yard  data  (fig.  3;  the  scaled  results  are  within  1  dB 

of  the  measured  values.  Scaling  with  ideal  acoustic  factors  gives 
*In  the  data  of  reference  5*  contributions  from  the  bubble  pulse 
were  predominant  for  frequencies  below  about  300  cps,  as  is  in¬ 
dicated  on  figure  3.  Bubble  pulses  were  not  considered  in  the 
digital  coaqjutatlons  represented  by  solid  lines  in  figures  3  and 

4. 
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spectra  that  are  too  high  at  the  higher  frequencies.  As  would 
be  expected  from  figures  1  and  2,  the  errors  are  greater  for  the 
larger  charge  and  longer  range  of  figure  4  than  for  the  scaled 
conditions  of  figure  3,  In  figure  4,  although  the  values  scaled 
with  ideal  acoustic  values  are  only  about  2  dB  too  high  at  1  kc, 
the  error  increases  to  about  +10  dB  at  a  frequency  of  10  kc. 

4.  ALLOWABLE  RANGE  OF  VARIABLES 

4.1  The  explosion  wave  scaling  factors  described  here  are 
applicable  only  so  long  as  the  shock  wave  exhibits  the  non¬ 
linear  characteristics  on  which  the  scaling  laws  are  based.  The 
ranges  and  weights  to  which  this  restricts  these  curves  are  the 
ranges  and  weights  for  which  the  explosion  similitude  equations 
describe  the  range  dependence  of  the  pressure  and  time  constant 
of  the  shock  wave.  Although  the  similitude  functions  are  not 
applicable  right  at  the  charge  Burface>  experiments  have  shown 
that  by  the  time  the  shock  wave  has  progressed  outward  to  only 
about  7  charge  radii ,  the  peak  pressure  and  time  constant  begin 
to  follow  the  similitude  relationships  (ref.  (1)).  The  real 
question  is  how  far  out  from  the  charge  these  same  relationships 
are  valid.  Arons  (ref.  (4))  found  that  the  l/ftl.13  pressure 
decay  and  the  time  constant  function  used  here  could  be  used 

to  describe  the  data  down  to  values  of  W1/?^  (cube  root  of 
weight  in  pounds  and  range  in  yards)  of  about  3.33  x  10-4.  Only 
this  one  set  of  data  is  available  for  such  small  values  of  wV3/fc 
and  the  experimental  scatter  Is  great.  Furthermore,  it  la  clear 
that  at  some  low  value  of  Vfl/3/fo  the  exponents  of  the  similitude 
equations  must  approach  the  acoustic  values;  Bhoek  wave  theory 
predicts  that  they  would  have  done  bo  at  shorter  ranges  than 
experiment  indicated.  In  view  of  these  considerations  it  is 
recommended  that  the  explosion  scaling  laws  developed  here  be 
restricted  to  the  interval  4  x  10“4  s  w W^/R  «  3  where  w  is 
charge  weight  in  pounds,  R  is  range  in  yards. 

4.2  In  addition  to  the  restriction  on  charge  weights  and 
ranges  noted  above,  these  scaling  laws  also  carry  a  frequency 
restriction.  The  scaling  functions  developed  here  assume  a 
steep  fronted  (l.e.,  zero  time  of  rise  to  the  maximum  pressure) 
exponential  wave  and  may  overestimate  the  spectrum  level  at 
high  frequencies.  Although  the  shook  wave  is  steep-fronted  near 
the  charge,  as  the  wave  propagates  outward  the  front  becomes 
increasingly  rounded  due  to  dissipation.  At  the  same  time, 
there  is  a  finite  amplitude  effect  which  tends  to  sharpen  up 

the  front  of  the  propagating  wave.  Arons  has  considered  the 
combined  results  of  these  effects  and  gives  theoretical  values 
of  rise  time  vs  range  for  various  charge  sizes  in  figure  13  of 
reference  (6) .  The  influence  of  the  rise  time  on  the  shock  wave 
spectrum  is  shown  for  several  hypothetical  curves  in  reference  (7). 
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From  references  (6)  and  (7)  it  appears  that,  for  the  stated 
Interval  of  Wv3/fa,  spectra  obtained  with  the  functions  shown 
here  are  probably  not  seriously  in  error  for  frequencies  below 
about  20  kc. 

4.3  The  lower  limit  of  frequencies  for  which  these  scaling 
laws  are  suitable  is  also  dependent  upon  the  exponential  approxi¬ 
mation,  which  neglects  the  slowly-decaying  pressure  of  the  shock 
wave  tail.  To  compensate  in  part  for  the  resulting  loss  of  low- 
frequency  components,  Weston  (ref,  (5))  used  an  average  value 

of  time  eonstant  (e  13  impulse/peak  pressure)  for  computing  his 
theoretical  spectrum  of  a  one-pound  charge  at  100  yards.  Such 
treatment  would  not  be  feasible  for  the  scaling  methods  shown 
here,  since  the  correction  to  the  similitude  e  that  would  be 
required  would  vary  depending  upon  the  charge  weights  and  ranges 
of  interest.  Furthermore,  for  charges  fired  at  fairly  shallow 
depths,  the  low-frequency  portion  of  the  spectrum  is  subject  to 
considerable  variation  due  to  bubble  pulses  and  surface  reflect¬ 
ions,  which  are  functions  of  charge  and  receiver  depths.  Con¬ 
sequently,  for  general  use  it  seems  more  appropriate  to  use  the 
similitude  values  of  time  constant  with  the  exponential  approxi¬ 
mation,  and  restrict  the  frequency  range  over  which  the  functions 
are  employed.  For  small  charges  the  Bpectrum  computed  for  the 
exponential  approximation  agrees  with  the  measured  spectrum  with¬ 
in  2  dB  or  better  for  frequencies  above  about  500  cps  (see 
Appendix  1).  Since  the  entire  pressure  time  curve,  as  well  as 
the  peak  pressure,  can  be  scaled  in  terms  of  wV3/tt,  one  would 
expect  a  similar  relationship  for  other  weights  and  ranges . 
Consequently,  500  cps  is  taken  as  a  reasonable  lower  frequency 
limit  for  the  scaling  relationships  given  here. 

4.4  To  summarize  the  above  restrictions,  the  scaling  laws 
given  here  should  not  be  used  outside  the  following  ranges  of 
variables : 

4  x  lO-*  ,  WV3/R  (lb  VS/Vd)  X  3 
0.5  <  frequency  (kc)  <  20 

Some  small  errors  due  to  the  assumptions  used  may  be  present  at 
the  extreme  values  of  these  ranges. 

5.  SUMMARY 

5.1  Scaling  factors  with  which  the  underwater  shock  wave 
energy  spectrum  for  one  set  of  experimental  conditions  (l.e., 
explosive  composition,  charge  weight,  and  range)  can  be  scaled 
to  another  set  of  conditions  are  derived  by  combining  shock  wave 
similitude  equations  with  the  Fourier  integral  for  the  spectrum 
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of  an  exponential  pulae.  These  factors  are  given  by  Equations 
(4)  and  (5). 

5.2  Scaling  factors  are  shown  graphically  for  the  case  of 
constant  range  with  varying  charge  weight  (fig.  1)  and  the  case 
of  constant  charge  weight  with  varying  range  (fig.  2).  Examples 
of  the  use  of  these  graphs  are  included. 

5.3  For  comparison  with  shock  wave  scaling  factors,  analogous 
scaling  factors  for  an  ideal  acoustic  wave  (i.e.,  one  in  which 
the  pressure  decay  with  distance  is  1/ft  and  the  wave  shape  is 
constant)  are  also  shown  in  figures  1  and  2.  For  charge  weights 
between  about  one  and  100  pounds,  and  for  ranges  between  about 

10  and  1000  yards,  there  is  little  difference  between  these  two 
sets  of  scaling  factors.  If  shorter  ranges  and  smaller  charge 
weights  are  involved,  however,  significant  errors  can  result  if 
shock  wave  spectra  are  scaled  with  the  ideal  acoustic  wave 
factors . 


5.4  Spectra  obtained  from  digital  analysis  of  shock  waves 
measured  at  one  yard  from  a  one-pound  charge  and  at  280  yardB 
from  a  1.8-pound  charge  are  shown  in  figures  3  and  4.  Data  given 
in  reference  (5)  for  octave  band  measurements  at  100-yards  from 

a  one-pound  charge  are  also  shown  in  figure  3.  The  measured 
one -pound, one -yard  spectrum  is  scaled  to  the  other  two  setB  of 
conditions  with  both  the  explosion  factors  and  the  ideal  acoustic 
factors  of  figures  1  and  2, 

5.5  Spectrum  scaling  with  explosion  wave  factors  gives 
better  agreement  with  the  data.  Values  scaled  with  ideal  acoustic 
factors  are  too  high  at  higher  frequencies,  with  worse  agreement 
for  the  longer  range  and  larger  weight  case  (fig.  4).  In  scaling 
from  the  one-pound,  one-yard  data  to  1.8-pound,  280-yard  levels, 
the  error  with  ideal  acoustic  scaling  is  only  about  +  2  dB  for 
frequencies  below  1  kc,but  increases  to  +  10  dB  at  a  frequency 

of  10  kc. 


5.6  The  allowable  range  of  variables  for  the  spectrum 
scaling  functions  given  here  is  estimated  to  be: 

0.5  <  frequency  (kc)  <  20 

4  x  10"^  s  £  3i  where  W  is  charge  weight  in  pounds 

and  R  is  range  in  yards. 

5.7  Spectral  analysis  of  the  measured  shock  waves  is  dis¬ 
cussed  in  Appendix  1. 

5.8  A  method  for  estimating  the  effect  of  charge  composition 
on  spectrum  level  is  given  in  Appendix  2, 
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CHARGE  WT (LB) 

(OR  RATIOS  OF  CHARGE  WEIGHTS) 

FIG  I  FACTORS  FOR  SCALING  ENERGY  SPECTRUM  FOR  I  "LB  TNT 
CHARGE  TO  OTHER  CHARGE  WEIGHTS,  SAME  RANGE 
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RANGE  FROM  CHARGE  (YARDS) 

(OR  I02X  RATIOS  OF  RANGES) 

FIG  2  FACTORS  FOR  SCALING  ENERGY  SPECTRUM  OF  GIVEN  WEIGHT 
TNT  CHARGE  FROM  IOO  YARD  TO  OTHER  RANGES 
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FIG.  3  ENERGY  FLUX  SPECTRA  FOR  I  LB  TNT  CHARGES 

AT  I  YARD  AND  IOO  YARDS 
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FROM  A  1.8  POUND  TNT  CHARGE 
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APPENDIX  1.  SPECTRAL  ANALYSIS  OP  MEASURED  SHOCK  WAVES 


1.  The  energy  spectra  of  several  underwater  pressure-time 
curves  recorded  at  various  distances  from  small  point  charges 
were  computed  by  Fourier  analysis,  using  an  IBM  7090  machine 
program  developed  for  use  with  arbitrary  wave  forms.  The 


functions  evaluated  were: 

a<*5  - } 

\ 

2 

N+l  .  K+l 

£  p  (T)  sin  os  T  +  £  p  (T)  cos  co  T 
n*l  1  jn**l 

2 

M  - 

E(f)  -|0  U(f) 


where  T  «  (n  -  |  and  <0  *>  2mf 

t  Is  the  total  time  of  integration  (seconds);  N  is  the  number  of 
equal  time  increments  in  the  Interval  t:  u>  Is  the  angular  frequen¬ 
cy  (radians)  and  p(T)  is  the  amplitude  (dyneB/cm2)  of  the  pressure- 
time  curve  at  the  mid-point  of  the  nth  time  increment. 


2,  In  the  energy  calculations  sn  average  value  of 

4  2 

pc  -  15.25  x  10  grams/cm  /sec 

was  used.  Variations  of  water  temperature  between  -2°C  and  28°c, 
salinities  between  34°/oo  and  $6°/00 ,  and  densities  between  1.0 
g/cm3  and  1.05  g/em3  result  in  a  spread  of  only  1.74  x  104 
g/cnr/sec  in  the  values  of  pc. 

3.  Spebtra  obtained  from  such  point-wise  integration,  which 
takes  into  account  the  slowly-decaying  tall  of  the  shock  wave, 
differ'jjrom  theoretical  spectra  for  the  initial  exponential 
portion  of  the  pulse  only  at  frequencies  below  about  1000  ops 
for  the  conditions  measured  here.  At  a  frequency  of  about  500 
cpa  the  spectrum  for  the  exponential  approximation  is  down  2  dB 
from  the  measured  spectrum;  the  difference  increases  as  frequency 
decreases . 


4.  The  spectrum  at  1  yard  from  a  1-pound  charge  (fig.  3  of 
report)  was  computed  from  an  oscillograph  recording  made  about 
one  yard  away  from  a  one-pound  charge  in  the  Potomac  River*; 
both  charge  and  tourmaline  gage  were  suspended  at  mid-depth  In 


Measurements  were  actually  made  at  2.5  feet  from  the  charge  and 
the  spectrum  was  then  scaled  to  a  one-yard  range  using  Equation 
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about  30  feet  of  water.  The  record  for  a  1.8-pound  charge  (fig. 
4  of  report)  was  obtained  from  an  oscillograph  playback  of  a 
tape  recording  made  at  sea.  In  this  case  the  charge  was  fired 
at  a  depth  of  356  feet  and  recorded  on  a  hydrophone  296  feet 
deep  and  844  feet  (slant  range)  from  the  charge.  The  frequency 
response  of  the  oscilloscope  recording  equipment  was  good  from 
essentially  zero  to  100  kc;  with  the  tape  recording  equipment 
the  frequency  range  was  only  30  ops  to  20  ko.  The  measured 
spectrum  for  the  latter  recor'd  (fig.  4)  reflects  the  poor  low- 
frequency  response  of  the  system. 

5.  For  both  of  these  shock  wave  records  integration  was  carried 
out  over  the  time  Interval  between  the  initial  pressure  dis¬ 
continuity  at  the  shock  wave  arrival  and  the  subsequent  decay 
of  the  wave  to  approximately  zero  pressure.  The  integration 
Intervals  actually  used  were  1.25  millisecs  for  the  one-pound 
charge  and  I.29  millisecs  for  the  1.8-pound  charge.  Sllfko 
(ref.  10)  has  shown  that  for  an  exponential  wave  the  interval 
of  integration  has  little  influence  on  the  spectrum,  provided 
the  Interval  is  at  least  as  long  as  five  times  the  decay  con¬ 
stant  of  the  wave.  Here  the  integrations  were  extended  over 
about  10  and  4  times  the  decay  constants  of  the  exponential 
portions  of  the  pulse  for  the  1-pound  and  1.8-pound  charges, 
respectively. 
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APPENDIX  2 

ESTIMATING  THE  EFFECT  OF  CHARGE  COMPOSITION  ON  SPECTRUM 
LEVEL  FOR  CONSTANT  CHARGE  WEIGHT  AND  RANGE 


1,  Although  TNT,  or  som«  explosive  material  that  la  very  similar 
in  performance  and  propagation  characteristics,  is  usually  used 
for  underwater  explosive  sound  sources,  the  effect  of  explosive 
composition  Is  of  Interest  when  new  explosive  materials  are 
considered  for  such  applications*  When  the  similitude  constants 
for  any  explosive  material  are  known  the  spectrum  scaling 
factors,  auoh  as  those  shown  for  TNT  In  figures  1  and  2  of  the 
report,  oan  be  determined  from  Equations  (5)  and  (5)  of  the 
report.  In  many  oases,  however,  one  is  primarily  Interested  In 
an  estimate  of  the  approximate  difference  between  the  spectrum 
of  the  new  material  and  that  of  a  known  material,  say  TNT,  for 
the  same  charge  weight  and  range.  For  such  comparisons  the 
following  equations  relating  e,  the  change  In  shock  wave  energy 
spectrum  level  (in  dB),  to  the  similitude  parameters  (Equations 
(2)  and  (3)  of  report)  for  the  two  materials  are  useful: 


c  -  10  log  E2(f)  -  10  log  E^f)  •  C  +  Xg  +  Xf  (A.l) 

where  C  -  20  log  (k^J/tk^)  (A  .2) 

3%  -  20  (a2  -  aj  +  82  "  *i)l°8  (V1^/*)  (A .3) 

Xf  -  10  log  [(1  +  a2  e^j/d  +  a2  §22)]  (A  .4) 

and  Wj  -  Wg,  R^  ■  Rg. 


The  quantity  C  is  Indicative  of  the  relative  energies  of  the 
two  materials  and  Is  Independent  of  charge  weight  and  range. 

The  variables  Xg  and  Xf  are  dependent  upon  the  similitude,  or 
propagation,  characteristics  of  the  shock  waves  of  the  two 
materials;  the  charge  weight  and  range  for  which  the  two 
materials  are  being  compared  enter  directly  Into  the  X~.  factor, 
and  In  the  9  terms  of  the  Xf  factor.  L 

2.  An  approximate  comparison  of  the  two  materials  can  usually 
be  made  very  easily  from  these  equations.  In  practice  It  has 
been  found  that  the  sum  of  the  similitude  exponents  (a2  -  ai 
+  go  "  Bi)  In  Equation  (A. 3)  is  very  nearly  zero  for  most  1 
hlgn  explosive  materials,  so  that  Xg  is  negligibly  small.  In 
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addition,  for  charge  weights  of  no  more  than  a  few  pounds,  Xf 
for  extreme  values  of  frequency  can  be  approximated  by: 

Xf  ■  0  for  frequencies  below  about  1  kc 

Xf  «  20  log  (e^/Bg)  for  frequencies  above  about  8  kc 

so  that  the  effects  of  this  term  on  spectrum  level  can  be 
estimated  quickly  for  several  values  of  wV3/ft. 


3.  For  the  occasional  oase  in  which  more  detailed  comparisons 
are  needed,  graphs  of  the  functions  Xg  and  Xf  are  convenient. 


Two  such  graphs,  covering  the  range  of  variables  customarily 
encountered  with  high  explosive  underwater  ohargeB,  are  shown 
in  figures  A.l  and  A. 2.  In  figure  A.l  the  quantity  Xg  is 
plotted  vs  wl/3/R  with 


r\  ■  (a2  -  a-^  +  Bg  “  B],)  (A. 5) 

as  parameter.  Note  that  in  this  graph  the  range  (R)  is  used 
in  feet,  as  is  the  custom  for  explosion  similitude  curves, 
rather  than  in  yards,  as  in  the  body  of  the  report.  In  figure 
A. 2,  Xf  vs  a)  e1  is  shown  for  several  values  of  Bg/Oi*  The  U8C 

of  Equation  (A.l)  and  figures  A.l  and  A. 2  for  ®°®P{|rln8 
energy  spectrum  for  some  hypothetical  explosive  Y  ,  which  is 
uncommonly  different  from  TNT,  with  that  of  TNT  is  illustrated 
below. 

4.  Assume  the  explosive  "Y",  for  which  the  subscript  "2"  is  used, 
has  the  following  similitude  values: 

p2  •  1.10  x  104  peak  pressure  in  psi  (A. 6) 

e2  -  0.10  («VV)-°-09  K1/3,  time  constant  in  milllsec(A .7) 

Here  W  is  in  pounds  and  R  in  feet.  Prom  reference  (4)  the 
corresponding  equations  for  TNT,  which  has  the  subscript  1  , 

are: 

Pa  -  2.16  x  104  (W1/3/^)1*13  (a .8) 

#1  .  0.058  (W1/3/*)’0,22  wV3  (A. 9) 
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The  comparison  will  be  made  for.  equal  ranges  and  weights  at  two 
values  of  W  viz.,  (1)  wl/3/R  •  1.0,  whloh  Is  quite  near 

the  charge  and  might  be  thought  of  as  an  approximate  comparison 
?L"°^rc!  lavel8»  «nd  (2>  VW/R  -  10-3,  which  corresponds  to 
1000  feet  from  a  one-pound  charge  and  Indicates  the  differences 
in  propagation  of  the  two  shock  waves.  The  spectrum  levels  will 
be  compared  at  two  frequencies,  1  and  8  kc. 


Prom  Equation  (A. 2), 


C  -  20  log  (kgtttgJ/Xkjn^) 


oft  1 1.10  x  0.10 

20  102  ?7T6'""x"  "6.658 


1.1  dB 


From  Equation  (A. 5), 

ss?  ‘-1*  *■ 


rj  *  1.05  - 
-  0  dB  at  wl, 


0.09  +  0.22  -  +  0.05i 

1.0  and  Xj  »  -3  dB  at 


For  computing  the  quantity  Xf  the  actual  oharge  weight  must 
be  used,  since  it  appears  in  theAequation  for  g.  If  a  oharge 
weight  of  one  pound  is  assumed  then: 

CW  m  l.o,  fig  «•  10  4  sec  and  0^  ■  0.58  x  10”4  sec; 

62/6!  «  1.725 

•  wV^/R  „  I0~3f  «  1.86  x  10“4  sec  and  ex  -  2.65  x  10'4  tec; 

62/^1  —  0.702 


and  with  figure  A. 2  the  Xj  values  tabulated  below  can  be  read. 


f  (kc)  <a 

1  0.63  X  104 

8  5.03  x  104 

Summarizing  the 
for  «  in  dB  are: 


wV3/h  ■  1.0  w l^Ai  -  io"3 


CJ0J 

Xf(dB) 

oej 

Xf(dB) 

0.365 

-1 

1.67 

+2 

2.92 

-4 

13.3 

+3 

above 

factors 

(Equation  A.l) 

the  value 

f  -  1  kc 
f  -  8  kc 


W Xft/R~  l.Q 


W^/ft  ■  io"3 


c  —  t2  dB  c  ■  -2  dB 

c  -  -5  dB  c  -  -1  dB 
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Thus,  very  near  the  charge  the  energy  spectrum  level  of  a  one- 
pound  charge  of  explosive  "Y"  1b  below  that  of  TNT,  the 
difference  increasing  from  -2  dB  at  1  kc  to  -5  dB  at  o  kc.  Be¬ 
cause  of  differences  in  the  propagation  laws  for  the  two 
materials,  however,  at  1000  feet  from  the  burst  the  relationship 
of  the  two  spectra  haa  changed;  here  the  level  for  explosive  Y 
is  still  2  dB  below  that  of  TNT  at  the  lower  frequency  (1  kc) 
but  only  1  dB  below  the  TNT  values  at  8  kc. 
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